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Abstract: We present an overview of possible imprints of non-extensitivity in particle and nucler physics. 
Special emphasis is put on the intrinsic fluctuations present in the system under consideration as the 
possible source of nonextensivity. The possible connection of nonextensivity and the self organized criticality 
apparently being observed in some cosmic rays and hadronic experiments will also be discussed. 
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1 Introduction - how we got to it 

Our encounter with notion of nonextensivity originated with ob- 
servation that in some cosmic ray data (like depth distribution 
of starting points of cascades in Pamir lead chamber |Tj) one en- 
counters deviations from the usual exponential distributions of 
some variables towards the power-like ones: 
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Here N denotes the number of counts at depth T (cf. [JTJ for de- 
tails). Whereas in QT]] we have proposed as explanation a possible 
fluctuations of the mean free path A in eq. ([]]) characterised by 

relative variance to = ^ — — - > 0.2, in H the same data were 
fitted by power-like (Levy type) formula as above keeping A fixed 
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and setting q = 1.3. In this way we have learned about Tsallis 
statistics and Tsallis nonextensive entropy and distributions [f3j. 
Following this approach we have then demonstrated [j||,|5j] that, 
indeed, 



with g=l±cjforg>l (+) and q < 1 (— ), i.e., there is con- 
nection between the measure of fluctuations u and the measure 
of nonextensivity q (it has been confirmed recently in [J6J). This 
will be the first subject discussed here. The other will be con- 
nected with the self-organized criticality form of nonextensivity 
as apparently seen in hadronic processes. 

2 Temperature fluctuations? 

The most interesting example of the possible existence of such 
fluctuations is the trace of power like behaviour of the trans- 
verse momentum distribution in multiparticle production pro- 
cesses encountered in heavy ion collisions |Q,|8j]. Such collisions 
are of special interest because they are the only place where new 
state of matter, the Quark Gluon Plasma, can be produced 
Transverse momentum distributions are believed to provide in- 
formation on the temperature T of reaction, which is given by 
the inverse slope of dN/dpx, if it is exponential one. If it is not, 
question arises what we are really measuring. One explanation is 
the possible flow of the matter, the other, which we shall follow 
here, is the nonextensivity (or rather fluctuations leading to it). 
Namely, as was discussed in detail in |7|] the extreme conditions of 
high density and temperature occuring in ultrarelativistic heavy 
ion collisions can invalidate the usual BG approach and lead to 
q > 1, i.e., to 
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Here m is the mass of produced particle and T is, for the q = 1 
case, the temperature of the hadronic system produced. Although 
very small (\q— 1| ~ 0.015) this deviation, if interpreted according 
to eq. (0)), leads to quite large relative fluctuations of tempera- 
ture existing in the nuclear collisions, AT/T ~ 0.12. It is impor- 
tant to stress that these are fluctuations existing in small parts 
of hadronic system in respect to the whole system rather than of 
the event-by-event type for which, AT/T = 0.06/ ^/N — > for 
large N (cf. [|nj] for relevant references). Such fluctuations are 



potentially very interesting because they provide direct measure 
of the total heat capacity C of the system: 

— ^ = c = u = q ~ (4) 



(f3 = ^) in terms of u = q — 1. Therefore, measuring both the 
temperature of reaction T and (via nonextensivity q 1) its total 
heat capacity C, one can not only check whether an approximate 
thermodynamics state is formed in a single collision but also what 
are its theromdynamical properties (especially in what concerns 
the existence and type of the possible phase transitions [pp . 

To observe such fluctuations an event-by- event analysis of data is 
needed [[TT| . Two scenarios must be checked: (a) T is constant in 



each event but because of different initial conditions it fluctuates 
from event to event and (b) T fluctuates in each event around 
some mean value Tq. Fig. 1 shows typical event obtained in sim- 
ulations performed for central Pb + Pb collisions taking place for 
beam energy equal Eb ea m = 3 ATeV in which density of particles 
in central region (defined by rapidity window —1.5 < y < 1.5) is 
equal to ^ = 6000 (this is the usual value given by commonly 
used event generators [|T2|). In case (a) in each event one expects 
exponential dependence with T = T event and possible departure 
from it would occur only after averaging over all events. It would 
reflect fluctuations originating from different initial conditions for 
each particular collision. This situation is illustrated in Fig. la 
where pr distributions for T = 200 MeV (black symbols) and 
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T = 250 MeV (open symbols) are presented. Such values of T 
correspond to typical uncertainties in T expected at LHC ac- 
celerator at CERN. Notice that both curves presented here are 
straight lines. In case (6) one should observe departure from the 
exponential behaviour already on the single event level and it 
should be fully given by q > 1. It reflects situation when, due to 
some intrinsically dynamical reasons, different parts of a given 
event can have different temperatures ||,|5[] . In Fig. lb black sym- 
bols represent exponential dependence obtained for T = 200 MeV 
(the same as in Fig. la), open symbols show the power-like depen- 
dence as given by (0) with the same T and with q = 1.05 (notice 
that the corresponding curve bends slightly upward here). In this 
typical event we have ~ 18000 secondaries, i.e., practically the 
maximal possible number. Notice that points with highest px cor- 
respond already to single particles. As one can see, experimental 
differentiation between these two scenarios will be very difficult, 
although not totally impossible. On the other hand, if successful 
it would be very rewarding - as we have stressed before. 

One should mention at this point that to the same cathegory 
of fluctuating temperature belongs also attempt |T3j] to fit en- 



ergy spectra in both the longitudinal and transverse momenta 
of particles produced in the e + e~ annihilation processes at high 
energies, novel nonextensive formulation of Hagedorn statistical 
model of hadronization process |TD|,|5f and description of single 



particle spectra |F| 



3 Nonexponential decays 



Another hint for intrinsic fluctuations operating in the physical 
system could be the known phenomenon of nonexponential de- 
cays [|T5j. Spontaneous decays of quantum-mechanical unstable 
systems cannot be described by the pure exponential law (nei- 
ther for short nor for long times) and survival time probability 
is P(t) oc t~ s instead of exponential one. It turns out |R3[] that 
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by using random matrix approach, such decays can emerge in a 
natural way from the possible fluctuations of parameter 7 = 1/r 
in the exponential distribution P(t) = exp(— 7^). Namely, in the 
case of multichannel decays (with v channels of equal widths 
involved) one gets fluctuating widths distributed according to 
gamma function 
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and strength of their fluctuations is given by relative variance 

— ^—^2 — = which decreases with increasing z/. According to 
, it means therefore that, indeed, 
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with the nonextensivity parameter equal to q = 1 + ^. 



4 Self-organized criticality in cascade processes 



For non-equilibrium phenomena the sources of power-law distri- 
butions are self-organized criticality (SOC - nonequilibrium sys- 
tems are continuously driven by their own internal dynamic to 
a critical state with power-laws omnipresent) [|TB| and stochastic 
multiplicative processes (power-law is generated by the presence 
of underlying replication of events) [PT7I1 - All they can be unified 
in terms of generalized nonextensive statistics 0] . It has been re- 
cently argued [JT^,[T^] that in systems of quarks and gluons formed 
in collision processes there exists evidence of SOC behaviour. The 
picture proposed in [|T8j is that system of interacting soft gluons 
can be, and should be, considered as an open, dynamical, com- 
plex system with many degrees of freedom, remaining in general 
far from thermal and/or chemical equilibrium. Applied first to 
description of formation of colour-singlet gluon clusters in inelas- 
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tic diffraction scattering this argumentation has been recently 
extended to cover also quarks This allowed to describe the 



existing data on high transverse momentum jet production in 
a uniform way (jet is, loosely speaking, a bunch of collimated 
hadrons going in one direction). Namely, the Et jet cross section 
(E T is the energy of such jet measured in direction perpendicular 
to the collision axis) is given by simple power law: 

S(T oc£?°, (7) 



dErdr] 



with a ~ 5, 7 and 9 for small, intermediate and large values 
of Et, respectively, indicating three different possible scenarios 
operating at different ranges of Et (according to [[H| ). 



Similar situation has been encountered in cosmic ray physics 
where energy spectra of particles from atmospheric cascades ex- 
hibit power-like behaviour, which can be described by Levy type 
distribution [PD|1 . For sufficiently large energy fraction xn allo- 



cated to the cascade with N generations it is given by 




P(x N ) oc ( 7:rT ] , a = — — -, (8) 



where q is nonextensivity index depending on the number of gen- 
erations N, q = ~ — fL ^— , (c is iV-independent parameter which 
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should be fitted, here c = 0.55). For events without cascading 
N = 1 and q = 1, i.e., one gets usual exponential distribution. 
For large N (in practice for N > 6) q — > 3/2 (and a — > 2 in (g)), 
which is limiting value available here for q (explaining why in 
such experiments one always observes a < 2). This result tells us 
that the exactly power-like behaviour of spectra is achieved only 
asymptotically, for long enough cascades. This has been clearly 
demonstrated in [^Dj on some experimental data. 



Coming back to eq. we argue that, as can be seen in Fig. 2, 
the same data can be described by a suitably modified equivalent 
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of eq. (|): 



d ' G |, = o = c^-[l-Erf(a)], (9) 



this time with only one exponent for all values of Et'- ol = 5.01 
and 4.9 for energies yfs = 546 and 1800 GeV considered in [|T9 



respectively. The correction factor to the power law comes from 
accounting for smearing out of the initial conditions of the cas- 
cade (not present in (g)) with a = (1 — «)(^) + -^hiAo and 
Aq = ypy (ratios of the actual energy of jet Et to its threshold 

energy E^). This smearing was assumed in the following log- 
normal form: P(A)d(ln A) = ^^exp — ( lnA 2 ^ A o) _ Parameter 5 
is dispersion of the smearing distribution and is set to be equal 
5 = 0.74 and 0.79 for the above respective energies. The value 
a ~ 5 in (j9]) emerges in a natural way from that of 2 in (|8|). The 
point is that the main interaction in hadronic collisions proceeds 
between quarks, which are produced in the quark-gluon QCD 
cascade process from the initial dressed valence quarks in nu- 
cleon. Therefore jets are produced by collision of two such quarks 
with energy fractions x and y and energy M of such collision is 
M = y/ z ■ s ~ yfz (where z = x-y and where s is initial invariant 
energy of reaction squared). Therefore distribution in z is of the 
form of convolution: P(z) = P(x) <S> P(y) ~ z~ 3 , which results in 

da 1 dcr c da , , c 

— = —^77 ~ M- 6 or — - - M~ 5 . 10 

dz MdM dM v ' 



Because Et ~ M eq . (|TU1) becomes eq.(|7|). Therefore value of a = 
5 is what really comes from the pure cascade. All deviations from 
it needed to fit data come from the initial conditions properly 
accounted for. 
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5 Summary 



There is steadily growing evidence that some peculiar features 
observed in particle and nuclear physics (including cosmic rays) 
can be most consistently explained in terms of the suitable ap- 
plications of nonextensive statistic of Tsallis. Here we were able 
to show only some selected examples, more can be found in [|5j]. 
However, there is also some resistance towards this idea, the best 
example of which is provided in [j2T[] . It is shown there that mean 
multiplicity of neutral mesons produced in p — p collisions as a 
function of their mass (in the range from = 0.55 GeV to 
Mr = 9.5 GeV) and the transverse mass tjit spectra of pions (in 
the range of tjit — l-i-15 GeV), both show a remarkable universal 
behaviour following over 10 orders of magnitude the same power 
law function C ■ x~ p (with x = m or x = uit) with P ~ 10.1 and 
P ~ 9.6, respectively. In this work such a form was just postulated 
whereas it emerges naturally in g-statistics with q = 1+1/P ~ 1.1 
(quite close to results of |T3j]). We regard it as new, strong imprint 
of nonextensivity present in multiparticle production processes. 
This interpretation is additionally supported by the fact that in 
both cases considered in |2Tj constant c is the same. Apparently 
there is no such phenomenon in AA collisions which has simple 
and interesting explanation: in nuclear collisions volume of in- 
teraction is much bigger what makes the heat capacity C also 
bigger. This in turn, cf. eq.(|]), makes q smaller. On should then, 
indeed, expect that qhadronic » Qnudear, as observed. 
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Fig. 1. (a) Normal exponential pt distributions (i.e., q = 1) for T = 200 MeV (black 
symbols) and T = 250 MeV open symbols), (b) Typical event from central Pb + Pb 
at Ebeam = 3 A-TeV (cf. text for other details) for = 200 MeV for q = 1 (black 
symbols) exponential dependence and q = 1.05 (open symbols). 
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Fig. 2. The high energy pp collisions data used in Ql9|] (cf. their Fig. 1) fitted by our 
eq. (pi). See text for details. 



10 



